Giant current fluctuations in an overheated single electron transistor 
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Interplay of cotunneling and single-electron tunneling in a thermally isolated single-electron tran- 
sistor (SET) leads to peculiar overheating effects. In particular, there is an interesting crossover 
interval where the competition between cotunneling and single-electron tunneling changes to the 
dominance of the latter. In this interval, the current exhibits anomalous sensitivity to the effective 
electron temperature of the transistor island and its fluctuations. We present a detailed study of 
the current and temperature fluctuations at this interesting point. The methods implemented allow 
for a complete characterization of the distribution of the fluctuating quantities, well beyond the 
Gaussian approximation. We reveal and explore the parameter range where, for sufficiently small 
transistor islands, the current fluctuations become gigantic. In this regime, the optimal value of the 
current, its expectation value, and its standard deviation differ from each other by parametrically 
large factors. This situation is unique for transport in nanostructures and for electron transport 
in general. The origin of this spectacular effect is the exponential sensitivity of the current to the 
fluctuating effective temperature. 

PACS numbers: 73.23.Hk,44.10.+i,72.70.+m 



I. INTRODUCTION 

By its statistical physics definition, the temperature 
of an open electron system in equilibrium is fixed, with 
a value equal to the temperature of the reservoir it is 
connected to. 1 The energy of the system can still fluc- 
tuate because of the constant exchange of energy with 
its surroundings. Out of equilibrium, the entire concept 
of a temperature becomes ill-defined. If the internal re- 
laxation in the system is strong enough, however, the 
electron energy distribution still follows the equilibrium 
Fermi function with some effective temperature that is 
determined by a balance of the energy currents flowing 
into and out of the system. 2 In this out-of-equilibrium 
situation the intrinsic fluctuations of the energy currents 
translate directly to fluctuations in the effective temper- 
ature of the system, provided that the electron-electron 
relaxation time is significantly shorter than the energy 
relaxation time to the reservoirs or the phonon bath. 3 If 
the electron-phonon interaction in the system is weak, 
the dominant energy flows are to the reservoirs which 
are used to drive the system out of equilibrium. In this 
case the system is said to be fully overheated. As the 
size of the system becomes ever smaller, full overheating 
becomes easier to achieve, or, harder to avoid. 

One example of a system where overheating is impor- 
tant to take into account is the single-electron transis- 
tor (SET). 4 When the size of the Coulomb island is de- 
creased ever further in the pursuit for sensitivity, the 
temperature of the island begins to be affected by the 
electronic heat flows to the surroundings, e.g., leads. 5 ~ 7 
In an earlier paper 8 we studied the fully overheated 
single-electron transistor, and found out that the inclu- 
sion of inelastic cotunneling, 9 along with the often dom- 
inant sequential tunneling of single electrons, is impor- 
tant and leads to peculiar overheating effects. The elec- 



tron transport through the overheated SET is divided 
into three regimes: cotunneling dominated, competition, 
and single-electron dominated regimes. We also found 
that the electric current noise in the crossover interval is 
mostly contributed by the frequency interval correspond- 
ing to typical time scale of temperature equilibration. 
This is because of the anomalously strong temperature 
sensitivity of the current in this interval, which makes 
the temperature fluctuations directly visible in the cur- 
rent (noise). 

In this Article we explore and quantify further the tem- 
perature fluctuations and the associated current fluctu- 
ations in an overheated SET. We study the SET un- 
der conditions of well-developed Coulomb blockade, so 
that the junction conductances Gt,l, Gt,r <C e 2 /h. 
We assume that the level spacing on the island is small 
compared to its temperature and the charging energy, 
5 1 <C T <C Ec = e 2 /2C, which is the case in metallic 
SETs. To have a clear model, we assume a vanishing 
temperature in the leads. 10 Since a possible asymmetry 
of the SET does not affect our results qualitatively, we as- 
sume Gt,l — Gt,r — Gt and define a dimensionless con- 
ductance gr = Grfi/e 2 <C 1. Our focus is mainly on the 
fluctuations around the crossover between the competi- 
tion and single-electron regimes, H ~Vc = (\/2— l)Vth, 
i.e., well below the zero temperature Coulomb blockade 
threshold, eV th = 2(E C - eV G ) (see Fig. 1). Note that 
in the actual experimental setup eVc must be replaced 
with a combination of the source, drain, and gate volt- 
ages together with the respective mutual capacitances. 11 
Here we lump all this into one gate voltage dependence. 
The electron-electron relaxation time, r e _ e , is assumed 
short compared to the characteristic time scale for the 
energy relaxation to the leads, te, which in turn is as- 
sumed small compared to the electron-phonon relaxation 
time, T e ^ p h- 12 These conditions ensure that the island is 
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FIG. 1. (color online) (a) A schematic diagram of a single- 
electron transistor, biased by voltage Vt- The charge on the 
island can be tuned with the gate voltage Vg- (b) Coulomb 
diamonds in a symmetric SET. Blue dashed line shows the 
threshold voltage Vc- Changing the gate voltage has the effect 
of changing Vc- 



fully overheated and in quasi-equilibrium. 

Within the Coulomb blockade, the temperature of the 
island is fixed by a balance between energy flows in and 
out of the island due to inelastic cotunneling and sequen- 
tial tunneling (see Fig. 2). Inelastic cotunneling leaves 
behind an electron-hole excitation on the island, always 
heating it up. Sequential tunneling can either cool or 
heat the island, depending on bias voltage. As discussed 
in Ref. 8, there are three qualitatively different regions 
of bias voltages, revealed by the behavior of the average 
temperature: At the lowest voltages, sequential tunneling 
is completely suppressed and the average temperature is 
fixed by the inelastic cotunneling. Close to the Coulomb 
blockade threshold (and above), inelastic cotunneling ef- 
fects can be disregarded and sequential tunneling sets 
the temperature. The most interesting region is the in- 
termediate regime where the competition between these 
two processes gives rise to both anomalously sensitive 
temperature dependence of the current and to strongly 
non-Gaussian temperature fluctuations. In this regime, 
we find (see Sec. II. A and App. A) for the average tem- 
perature (T) w Tc where Tc satisfies 
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FIG. 2. (color online) Schematic energy balance of the single- 
electron transistor within the Coulomb blockade regime, 
where the allowed charge state energies are outside the bias 
window. In this case the temperature of the island is fixed 
by the balance between sequential tunneling (solid red ar- 
rows) and inelastic cotunneling (dashed blue arrows), where 
the strength of the previous depends on the exponentially 
small (oc exp(— 5i5 + /fcsT)) amount of empty (filled) single- 
electron states below (above) the Coulomb blockade thresh- 
old corresponding to the addition (removal) of an electron to 
(from) the island. This leads to the extreme temperature sen- 
sitivity of the current, which is the main topic of this paper. 

a being a dimensionless coefficient, a w 0.1. For exam- 
ple, for c/t = 10~ 3 , Tc /Vc ~ 0.1 and for smaller values 
of <3>T) Tc /Vc ~ l/[V21n(l/g r r)]. As we show below, the 
small value of Tc/Vc allows us to rigorously describe the 
anomalously large and strongly non-Gaussian tempera- 
ture and current fluctuations in an overheated SET. 

We find that the probability distributions of the tem- 
perature and current fluctuations are strongly non- 
Gaussian, and that the size of the transistor island has 
a large influence on these distributions. In particular, 
we find that for small islands with relatively large single- 
electron level spacing, 

'■*{%)'«>■ 

the mode, expectation value, and standard deviation 
of the temperature and current differ from each other 
by parametrically large factors (see Figs. 6 and 8, and 
Eqs. (32), (34) and (35)). This we show by first cast- 
ing the Keldysh action technique used in Refs. 3 and 8 
into the form of a Fokker-Planck equation, allowing us to 
describe the full normalized probability distribution, in- 
stead of the non-normalized large deviation function In P. 
For some realistic values, T c /Vc ~ 0.1 and Eq ~ lmeV, 
5i > 10 -8 eV, corresponding to roughly 10 9 atoms. 13 

This Article is structured as follows: Section II intro- 
duces the theoretical methods utilized in this paper. In 
Sec. Ill we concentrate on the regime of small fluctuations 
and in Sec. IV on fluctuations in systems with relatively 
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large single-electron level spacing. We conclude and dis- 
cuss our findings in Sec. V. In most of the text, we use 
units such that h = ks = e = 1, but restore them in the 
important results. 



II. APPROACH 

We use the extended Keldysh action technique 14 as the 
basis of our theoretical calculations. Within this frame- 
work it is straightforward to calculate the full statistics 
of, for example, electric and energy current. 15 ' 16 The ef- 
fective action of the system is augmented with counting 
fields XL(R)(t) and ^L(R)(t) for charge and energy transfer 
to the left (right) reservoir, respectively. The probabil- 
ity to obtain a certain realization of fluctuating energy 
currents is 



P[H L (t)]P[H R (t)] = / Va(tM R (t) 



x exp j- J dt (i£, L H L + \$ R H R - S(Z l ,Zr, £)) | , 

(2) 

where <S is the effective action describing the transport 
processes. 17 Analogous expression gives the probabilities 
for charge currents with the substitution £ i— > \. To 
ensure the conservation of energy, E, on the island, i.e., 
H L (t) + H R (t) — E(t), we use a delta functional, 18 

5[H L (t)+H R (t)-E(t)] = J T>£(t) 

x exp jiy* + (3) 

to generate the conditional probability 

P[H L (t),H R (t)} = J V^ L (t)V^ R (t)V^(t)VE(t) 

xcxp dt (i(£ L -Z)H L +i{£ R -Z)H R 

+^E-S{^ R ,E))\. (4) 



In the case of the single-electron transistor adding a delta 
functional for the conservation of charge is not necessary, 
since the effective action we use already conserves charge. 
Integrating over H L and H R yields the energy- conserving 
Keldysh partition function as 

Z = J Vt(t)VE(t)exp{- J dt (t(t)E(t)-S(t,E))y 

(5) 

where the imaginary unit has been absorbed into £. 

To calculate the probability of the island having energy 
between E* ...E* + AE during time r . . . r + At the 
functional integration over E(t) in Eq. (5) is performed 
over paths which satisfy this requirement. In general, the 



functional integration cannot be performed analytically. 
The variation of the exponent with respect to E and £ 
yields the semiclassical equations of motion, 



e = d^s, i = -a E s, 



(6) 



which, when solved with the desired boundary conditions 
E(t) = E* and £(— oo) — > 0, give the saddle point trajec- 
tories of E and £. The probability is then obtained as the 
saddle point value of the functional integral. 3 The sad- 
dle point approximation is valid when the fluctuations 
around the semiclassical trajectories are small. In our 
case this translates to the condition 5i <C (Tc /Vc) Vc- 
For another approach, we expand S to second order in 
£. Then the Gaussian integral over £ in Eq. (5) can be 
performed and the remaining functional integral trans- 
formed into a Fokker-Planck equation for the time evo- 
lution of the probability distribution of E: 14 



&P(E,t) 
dt 



= HV(E,t), 



d_ 
dE 



D 1 (E)V(E,t)-D 2 (E) 



dV(E,t) 



dE 



(7) 



Here, the operator H is obtained from <S(£, E) with the 
substitution £ —d/dE. Before applying this "quanti- 
zation rule", the operator H must be normally ordered, 
i.e., £ must be to the left of all E (see Appendix B). 
The Fokker-Planck equation also takes into account the 
Gaussian fluctuations around the saddle point trajectory, 
and is therefore more accurate than the saddle point ap- 
proximation. 

The stationary probability distribution satisfies 
HVst(E) = 0, and has the form 



V st (E) = const, x exp / dE 



ME) 
D 2 (E) 



(8) 



provided that / dEiD^E) / D 2 {E)) -> -oo when E -> 
±oo. The normalization condition J dE'P st (E) = 1 fixes 
the value of the constant prefactor. 

The time dependence of E can also be described in 
terms of a Langevin equation. The Fokker-Planck equa- 
tion is equivalent to a Langevin equation 19 



dE(t) 
dt 



Dim)) + 



dD 2 (E(t)) 
dE(t) 



V(t), (9) 



where rj(t) is a random "force" with mean (77(f)) = and 
variance (t](t)r](t')) = 2D 2 (E{t))S(t - t'). 

To obtain a probability distribution for the effective 
electron temperature, we must apply a model which de- 
scribes the relation between E and T. We first start by 
pointing out that the energy fluctuations show up as fluc- 
tuations of the electron distribution function /(e), where 
e is the energy of a given excitation. In a free-electron 
model, the total thermal energy on the island is related 
to this via 



E(t) 



J dee(/(e,f) 



/o(e)), 



(10) 
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where Si is the single particle level spacing and /o(e) is 
the zero-temperature electron distribution function. In 
the quasiequilibrium limit, /(e) tends to a Fermi distri- 
bution function, f(e,t) = (exp{e/[T(i)]} + with a 
fluctuating electron temperature. In this case we can in- 
tegrate Eq. (10) and obtain the usual expression for the 
free-electron heat capacity C(T), 



E : 



C(T)T 

2~ 



T 2 T 2 



6(5/ 



(11) 



Note that in principle the distribution of the fluctuating 
electron temperature may depend on the charge state n. 
However, under Coulomb blockade the relaxation rate of 
the excited charge states is larger than the temperature 
equilibration rate by a large factor, of the order of Eq/Si, 
and this effect can be disregarded in the following. 



A. Single-electron transistor 

For the overheated single-electron transistor the SE 
part of the action in the limit T -C Vb,Vc reads 8 



S sc =-g T Te- w ( T - 1+ V{ 



2 _ e (w+v b )i 

+ e w *}, (12) 



where W — (Vth — H)/2. In this article we concentrate 
on the cases where along the semiclassical trajectory, de- 
scribed by Eq. (6), £ <C l/V c , l/T c . In this case, a 
straightforward expansion to second order in £ yields 



S sc =-g T TV t l 1 e- w ' T 



W + 2TW + TV b - 



v; 



w 3 



TW Z - 



TV 2 



(13) 



Near the critical bias voltage, VJ, « Vc = (\/2 — l)14h, it 
is convenient to work with dimensionless variables 

t = V2T/V C , v = (V b - Vc)/V c , x = eV c , 
so that the action can be expanded to 

x 2 \ V c 



_ l-v/y/2 

-g T te * \(t-v)x- 



(14) 



valid when t <C 1, \v\ <C 1, \x\ <C 1. From this form we 
can immediately see that the relevant x w 6(t -ti)cl, 
justifying our expansion above. 

The cotunneling part of the action can be written as 8 



Scot —H cot £, + ^if^ot^, 



where H cnt and Si 



2 ^ 12 ' (15) 
i co t anu cot are the energy current and its fluc- 
tuations due to cotunneling, and a and /3 numerical fac- 
tors of the order of 0.1. Here we disregard the weak volt- 
age and temperature dependence of these terms. It turns 
out that the energy current noise due to cotunneling, i.e., 
the term proportional to /?, can be neglected: The single- 
electron term is proportional to gTtexp(— (1 — v/\/2)/t), 
whereas the cotunneling term is proportional to g 2 -,. Near 
the crossover the exponential term is of the order of grt~ 3 
(see Eq. (1)) and the single-electron contribution is there- 
fore g\t~ 2 > g\. 

The total action, including single-electron and cotun- 
neling contributions with the aforementioned approxima- 
tions, is then 

c + _ i = v/V5 ( x 2 \V C[ g\xV c 

S = -g T te * ( ( *~^ X ~~6~)t 2~ ' 



Using Eq. (11) we substitute 

d 6<5, d 



8E irHVgdt' 



(16) 
(17) 

6Sj ' \~ 125, ' ' " (18) 
the Fokker-Planck equation 

3^^(^) = -|{A(^(*.r)-D a (0|p(*,r)^ 



and obtain for the distribution of temperature, 

P(t,T) 



c -V [ E - c 



(19) 



where 



_ Wvg , ag T 8i l-v/y/2 
Di(t) — —(t — v)e * H : — ^ e * 



* 7T 2 Vc t 3 



D2(t) = 



5 1 1 l—p/%/2 
7 e ' ! 



ir 2 V c t 

which, using Eq. (8), yield the stationary distribution 



V st (t) oc exp ■ 



ir 2 V c 
Si 



=:t 3 



-t v — agrtexp 



V\/2 



+ ag T (l - u/\/2)Ei 



1 - v 



+ 



1 - v/V2 



(20) 
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Here, Ei(z) is the exponential integral function. 20 

To focus on the crossover regime, we introduce a new 
set of dimensionless parameters by changing the variables 
to 

t = V2T c /V c + 29(T c /Vc) 2 , 
v = V2T C /V C + 2u{T c /V c f, 

with the temperature at the crossover, Tc, defined 
through Eq. (1), see also Appendix A. We also intro- 
duce a dimensionless level spacing, 



V 



1 Si ( eV c 
W^eVc \k B T c 



characterizing the relative strength of the fluctuations, 
and a typical relaxation time for these fluctuations in 
the crossover interval, 



TV = 



2y/2n 2 h fk B T ( 



Sag 2 , (5/ \ eV c 



For a typical g T = 10~ 3 , r r w 1.7 x l(T 10 /<5/ eVs. When 
v, 9 <C Vc/Tc, and Tc/Vc <C 1, we can expand in the 
small parameter, Tc/Vc, and simplify the Fokker-Planck 
equation for the overheated SET considerably: 



d_ 



[-(0-v)e 6 + l-Ve e ]V{e,T) 



The stationary distribution becomes 
7> st (0)cxexp j-i (hp-v6 + e- 1 



(21) 



(22) 



which is the main result of this paper. The distribu- 
tion consists of two parts: the first, proportional to 1/T>, 
which is also present in the semiclassical (saddle-point) 
limit, and the second, independent of T>, describing the 
Gaussian fluctuations around the saddle-point trajecto- 
ries, which become significant as Si (and therefore T>) 
grows. Low-temperature fluctuations are exponentially 
suppressed since the cooling due to singe-electron tunnel- 
ing is also exponentially small. High-temperature fluctu- 
ations on the other hand do not have such suppressing 
mechanisms, hence they follow the usual Gaussian form. 



On the single-electron side of the crossover, but yet not 
far from the crossover point, 1 3> v ^> Tc/Vc, we can in 
addition neglect the exponential terms, proportional to 
g T , in Eq. (20): 



V st (t) oc exp 



Si 



1 s 

-t 3 



1 2 

-t 2 v 



(23) 



The distribution has a maximum at t = v. The cubic 
term favors low-temperature fluctuations, similar to that 
obtained in Ref. 3 for the noninteracting voltage biased 
island. 21 This is because far above the crossover SE pro- 
cesses serve as an efficient cooling mechanism. Similar 
to the results in Ref. 3, deviations from Gaussian statis- 
tics appear for temperature fluctuations of the order of 
average temperature, but their probability is greatly en- 
hanced since In V ~ T 3 /<5/Vj < T/Sj. 

In the competition regime, v < 0, we can approxi- 
mate the exponential integral with its asymptotic form, 
Ei(z) w e z (l/z + l/z 2 ), for large z. The resulting distri- 
bution is 



V st (t) oc exp 



7T 2 V C 

Si 

+ag T — 



1 

3 1 

e 



1 



t 2 v 



exp 



1 - v/V2 S 



(24) 



In contrast to the previous case, low-temperature fluctu- 
ations are exponentially suppressed for reasons explained 
at the end of Sec. II. 

Around the crossover, v < Tc/Vc, we have from 
Eq. (22) 



V st {9) oc exp <^ - 



1 fl 



V V 2 



v8 + e-° 



(25) 



This distribution is "half-Gaussian", i.e., Gaussian for 
high-temperature fluctuations but exponentially sup- 
pressed for low-temperature fluctuations. In this case 
the non-Gaussian character appears already for devia- 
tions of temperature of the order of ST ~ {Tc/Vc)Tq 
and their probability is even more enhanced, ln"P ~ 
{Tc/Vc)\T c /5i). 

The logarithm of these probability distributions is plot- 
ted in Figs. 3 and 4, which clearly show the half-Gaussian 
characteristics for voltages close to and below Vc- 



IV. GIANT FLUCTUATIONS 



III. REGIME OF SMALL FLUCTUATIONS 

In the limit of a metallic island, Pel, the last terms 
in the exponents in Eqs. (20) and (22) can be disregarded. 
In this limit the typical fluctuations are Gaussian and 
small: the variance is proportional to Si/Vc- It should 
be stressed that our approach still goes beyond the Gaus- 
sian approximation - we can also look at atypical, large 
fluctuations, which are non-Gaussian in character. 



In the previous section, we assumed that the fluctua- 
tions are small, P< 1. If V ~ 1, the typical fluctuations 
become non-Gaussian. For temperature, the fluctuations 
are still small, occurring at the scale of ST ~ Tc (Tc/Vc) ■ 
Since current is anomalously sensitive to temperature, 
these fluctuations manifest in much stronger fluctuations 
of the current. Throughout this section we work with the 
scaled variables, v and 9, and the results presented are 
valid for v, < V c /T c . 



G 




FIG. 3. (color online) Logarithm of the temperature fluctu- 
ation probability for some values of bias voltage v — (14 — 
V c )/Vc- In this plot a = 0.1 and g T = 10" 3 . 



FIG. 5. (color online) Probability distribution for the tem- 
perature of the island. Note how the maximum of the dis- 
tribution is shifted and the non-Gaussian tail becomes more 
prominent as T> is increased. 




FIG. 4. (color online) Logarithm of the temperature fluctu- 
ation probability for bias voltages v = —5 ... 5, bottom to 
top. 



A. Temperature fluctuations 
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FIG. 6. (color online) Expectation value of the temperature 
as a function of bias voltage. Dashed lines show the mode 
of 8, i.e., the maximum of the distribution. The two tend to 
each other for \u\ > T>. 



For a finite but small level spacing Sj, the distribution 
of temperature is given by Eq. (22) and plotted in Fig. 5. 
The maximum of the distribution is located at 

^={_4_~%), :<!! (26) 

The maximum is shifted to lower temperatures with in- 
creasing T> and it differs from the expectation value of 
the temperature when \v\/T> is small as shown in Fig. 6. 
The variance of 9 is plotted in Fig. 7. 

Analytical approximations for the expectation value 
and variance can be obtained in the limit \u\ ^> T>, v > 
and v < 0. In the first case the weight of the distribution 



is shifted to large 9, and we can neglect the term e in 
Eq. (22). In the latter case the weight is at small 9, and 
we can neglect the quadratic term in Eq. (22). We obtain 

V »C) = { -vU I < I • M » v - < 28 > 

For large positive or negative u, (9) ~ # max . Variance is 
proportional to T> as in the regime of small fluctuations. 
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FIG. 7. (color online) Variance of the temperature as a func- 
tion of bias voltage. For v/V 1 the variance approaches 
V. 



FIG. 8. (color online) Expectation value of the current as a 
function of bias voltage. Dashed lines show the mode of j, 
i.e., the maximum of the distribution. The ratio of these tend 
to exp(32?/2) for v/V > 1, V > 1. 



B. Induced current fluctuations 



In the overheated SET the electric current near the 
critical bias is given by 8 



1(6) = 



e 2 ag\e 2 Vc 
h 4k%T% ' 



(29) 



A fluctuating temperature leads directly to a fluctuating 
electric current. Since the current depends exponentially 
on the temperature, small fluctuations of temperature 
lead to large fluctuations of current. The Fokker-Planck 
equation for temperature is easily converted to a corre- 
sponding equation for the dimensionless current, j = e , 
with the prescription 

4=^, V(j,t) = -V(e = \nj,t), 
98 dj j 



resulting in 



(30) 



Note that this equation only includes the dominating con- 
tribution to the current fluctuations due to the tempera- 
ture fluctuations, and the intrinsic fluctuations (thermal 
and shot noise) arc disregarded. The stationary distribu- 
tion of current is then 



V s t(j) cx exp 



1 

V 



(Inj) 2 - fin j + 



21nj 



(31) 

Similarly to the case of temperature, the maximum of 
the distribution is shifted to lower values of j when D is 



increased, 
at 



The maximum of the distribution is located 



exp [v - 2D) . 
[-U + 2DY 1 



v>0 
u<0 



> 2D. 



(32) 



Due to the log normal character of the current distribu- 
tion, the expectation value and the most probable value 
of the current deviate as v is increased. This is shown in 
Fig. 8. For a small D, this deviation is proportional to D. 
The variance of j is plotted as a function of the (reduced) 
bias voltage in Fig. 9. Note that Var(j) is the variance 
of the instantaneous electric current, which is different 
from the zero frequency spectral noise power, more often 
encountered in the literature. 

To comprehend the unusual properties of the distribu- 
tion, let us first turn to the limit I)> 1, IH/f — 1- ^ n 
this case, the distribution reduces to a power-law func- 
tion 



^st(j) 



J 



1 

2-v/V ' 



(33) 



This form cannot be valid at all values of the current: 
It is cut off at small currents, j s ~ 1/22?, and crosses 
over to the log-normal distribution at large currents, % ~ 
exp(2?), j s C 1 € jb- Analytical approximations for the 
expectation value and variance can be obtained in the 
limit \v\ 3> D, v > and v < 0. Similarly to the case of 
the temperature distribution, in the first case the weight 
of the distribution is shifted to large j, and we can neglect 
the term 1/j in Eq. (31). In the latter case we can neglect 
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FIG. 9. (color online) Variance of the current as a function 
of bias voltage. Dashed lines show the square of the average 
current, {j) 2 /T>, for comparison. The ratio of these tend to 
exp(D) for v/V » 1, T> > 1. 
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FIG. 10. (color online) Three simulated realizations of current 
in an overheated SET at v = 0, T> = 3.3. The average current 
(for V -> 0) is j « 1.7. 



the (lnj) 2 term. We obtain 



(.?) = 
Var(j) = 



\v\ » V, 



exp(v-V/2), v>0 
-l/v, v<Q 

exp(2i/ - D)(exp(r>) - 1), v > 



(34) 

\v\ >X5. 
(35) 



We notice that for 1/ > 0, (j) parametrically exceeds 
ij) Umax = cxp(3X>/2), and ^J\ai{j) parametri- 
cally exceeds (j), ^/Var(j)/(j) = exp(X>/2). This signals 
a highly unusual distribution not satisfying, for instance, 
the central limit theorem conditions. 

We stress that the exponent of the power-law distri- 
bution can be readily tuned with bias voltage, dT4 = 
2Vc{Tc /Vc)&v ■ We are not aware of any other physical 
systems where a quantity exhibits a power-law distribu- 
tion with a tunable exponent. 

Using Eq. (9) allows us to write down a Langevin equa- 
tion for the time dependence of the current 



^j(i)+7(i(«))j (*)=#), 



(36) 



where = —j{v — \sxj + T>) — 1 and rj(t) satisfies 

(r?(t)) = and (r](t)r)(t')) = 2Vj(t) 3 5(t-t'). The current 
exhibits huge peaks at a low frequency that are orders of 
magnitude larger than the average current, as can be 
seen from Fig. 10. The time scale for these fluctuations, 
T r , is of the order of milliseconds for gr — 10~ 3 and 
<5/ = 0.01 K. 

Another way to access the time dependence of the cur- 
rent fluctuations is to study the saddle point equations 



for small fluctuations (6), i.e., vanishing 5], 

■ - 9 Q 

T T X- ^b, 

r r 6 = —S, 
ox 



where the action, 



'- v - 1 x- 



T 2 V 2 
> X V C 

12T 2 ' 



(37) 



(38) 



has been written in terms of the dimensionless, scaled, 
variables, v, 9, and x. As explained in Ref. 3, the 
saddle-point solutions follow the trajectory given by 
S(9(t),x s (6(t))) = with x s {8) ^ until a "measure- 
ment" at time r forces them to a trajectory with x$ = 0. 
As shown in Fig. 11, the shape of a fluctuation obtained 
from the saddle-point equations agrees with the Langevin 
result. This is expected, since for small T> the fluctua- 
tions around the saddle point trajectory are negligibly 
small. 



V. CONCLUSIONS 

We have studied the temperature fluctuations and the 
associated current fluctuations in an overheated SET. Fo- 
cusing on the crossover region between the competition of 
single-electron and cotunneling, and pure sequential tun- 
neling, we have found several interesting features, not 
seen in other commonly considered nanostructures. In 
overheated SETs with small islands the expectation value 
and the most probable value of the temperature differ 
from each other due to a tail in the probability distribu- 
tion of temperature, extending to significantly high tem- 
peratures. Due to exponential temperature sensitivity of 
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FIG. 11. (color online) Comparison of a saddle point fluctu- 
ation to the Langevin result at v — 0, T> — 0.02. 



the electric current, huge peaks, occurring at frequencies 
of the order of a few kilohertz, should be visible in an 
instantaneous measurement of the electric current. 

Experimentally, the challenge is to manufacture a tran- 
sistor island where, despite high resistance tunnel con- 
tacts, the electron-phonon heat current is negligible com- 
pared to the heat current to the leads. In earlier exper- 
iments on the measurement of noise in SETs (see for 
example Ref. 22), this has not been the case. According 
to our earlier evaluation 8 a transistor island with a vol- 
ume of V = 10 -4 /im 3 connected with gx = 10~ 3 tunnel 
contacts should allow one to detect these large current 
surges experimentally. 



Appendix A: Details on the rescaling 

The total energy current in the SET can be found from 
the action, Eq. (16), via 



H 



dS 

dx 



x=0 



_ 1-V/V2 Vc 
-g T te * [t-v)~- 



9t v c 



(Al) 



The heat balance follows from the requirement for the 
total energy current to vanish, leading to 

_ i—v/VS 

— te * (t — v) + agr = 0, 

v = -ag T t L e * + t. (A2) 

Next we make a change of variables to t = tc + 9t c and 
v = tc + vt C i where tc is for now an arbitrary number. 
If 9, v <^ 1/tc, we can expand in these small parameters 
to obtain 



Q _J ± 

-agrt c e*c -si +i 



(A3) 



Selecting 



ag^ c exp 



1 

71 



= l^t c 



1 



ln(l/ 5T ) 



< 1, 



simplifies the heat balance to 



(A4) 



(A5) 



and also implies tc < 1, thus justifying the expansion 
above for a large range of v and 9. We define the crossover 
point from the competition between single-electron and 
cotunneling to pure single-electron tunneling to be at 9 — 
0, i.e., t — tc = V^Tc/Vc, defining the temperature at 
the crossover. This happens when v = — 1, i.e., at a 
voltage of v = tc — t c . This rescaling simplifies also the 
Fokker-Planck equation, as seen in the main text, since 
it allows us to approximate 
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Appendix B: Relation between the Fokker— Planck 
equation and the functional integral 



Let us write the FP equation in the form 



d_ 
di 

d_ 
di 



/(«,*) 

fW,t)= I dqf(q,t) 



d 0^ 
- [£>i(«,t)/(3,*)] + ^2 Pa(9, *)/(«>*)] . 



-Dx (q, t) -5{q'-q) + D 2 (q, t)^8{q' - q) 



(Bl) 
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o o2 

S(q' -q)- AtD^q, t)-S(q' - q) + AtD 2 (q,t)-^8(q' - q) + 0(At 2 ) 



and integrate over time from t to t + At: 

f(q',t + At) = J dqf(q,t) 
Now we use the identity 5(q' — q) = J dp/2n exp(ip(q' — q)) to write 

/(</, t + At) = J dq^f(q, t)e^'-i) [l + ipAtD^q, t) - P 2 AtD 2 (q, t) + 0(At 2 )] . 

Iterating iV times so that (tjv — to)/ At = N gives 

N-i , ^ fN-i \ q q 

f(q N ,t N )= ]J / d( H^f(<lo,to)exp< ^2 At 



i=0 



At 



+ ip i D 1 (q i ,t i ) -p i D 2 (qi,U) 



which admits a functional integral representation 



f(q, t) = J VqVpf(q , t ) cxp dt [ipq + ipD^q, t) - p 2 D 2 (q, t)] J . 



(B2) 



(B3) 



(B4) 



(B5) 



This form of the functional integral is not unique (unlike Eq. (B4)), but depends on the discretization procedure. 23 
With the substitution ip — £, q E, we get a functional integral in the form of Eq. (5). This proves that the FP 
equation (7) corresponds to the partition function of Eq. (5). 
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